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Abstract 



The structure of time-dependent Gaussian solutions for the Kostin 
equation in dissipative quantum mechanics is analyzed. Expanding 
the generic external potential near the center of mass of the wave 
packet, one conclude that: the center of mass follows the dynam- 
ics of a classical particle under the external potential and a damping 
proportional to the velocity; the width of the wave packet satisfy a 
non-conservative Pinney equation. An appropriate perturbation the- 
ory is developed for the free particle case, solving the long standing 
problem of finding analytic expressions for square integrable solutions 
of the free Kostin equation. The associated Wigner function is also 
studied. 

Keywords: Kostin equation; dissipative quantum mechanics; damped 
Pinney equation. 

PACS: 02.30.Hq, 03.65.Ta, 42.50.Lc. 

1 Introduction 

Presently there is no universal consensus on how to quantize dissipative clas- 
sical systems, a subject of great practical and fundamental interest. Among 
the available alternatives, the Kostin logarithmic nonlinear Schrodinger equa- 
tion [9] provides a suitable modeling for dissipation in quantum mechanics. 
We note that presently there are no experiments measuring deviations, if 
any, from linearity in quantum theory. Hence, in spite of its drawbacks 
(breaking of the superposition principle, wrong modified frequency in the 
case of damped harmonic oscillations [6]), it is of interest to have a deeper 
understanding of the solutions of the Kostin equation. Moreover, recently 
the Kostin equation (also termed Schrodinger-Langevin equation) is attract- 
ing attention, in the context of dissipative time-dependent density functional 
theory [2U]- Similarly, the Kostin logarithmic term was suggested as an ap- 
propriate bath functional in time-dependent density functional theory for 
open quantum systems with unitary propagation [23] . In addition, some re- 
cent studies include a rigorous mathematical analysis of the Kostin equation 
in bounded domains [IU] , the stability analysis of the equilibrium states |22j , 
the search for numerical solutions under diverse confinement potentials [TT] . 
as well as connections on quantum combinatorial optimization, or quantum 
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annealing [1] . The existence of soliton-like solutions for Schrodinger-Langevin 
equations in a fluid form from the perspective of Nelson's stochastic mechan- 
ics [2] has been also investigated [13J. 

The interest on the Kostin model arises because: it preserves the norm 
of the wave function, in spite of the intrinsic nonlinearity; when written in 
the hydrodynamic form, one find a linear in the fluid velocity friction term, 
with a simple interpretation in terms of Ohmic dissipation; the Ehrenfest 
relations for dissipative systems are strictly satisfied. Further possible diffi- 
culties, relying on the existence of undamped stationary states, actually can 
be understood since these states might be interpreted as one- dimensional 
projections of the normal modes of motion of a system with many degrees of 
freedom coupled by conservative forces [7]. 

It is a long standing problem, to accurately describe the free particle 
dynamics in the Kostin model and so far only numerical solutions are avail- 
able [5], [7J. It is one of the main purposes of this work, to circumvent 
this difficulty by means of an appropriate perturbation theory in a damping 
dominated scenario. Furthermore, it was shown [2] that in the damped free 
particle case the Kostin equation can have no solitary wave type solutions of 
the form \1/ = + ct) case which lie in L 2 , where c is a constant. Also, 
in his original paper Kostin [H] have constructed some solutions of the plane 
wave like form \1/ = exp(i9(t)) exp[(i/h)(az(t)x + (3(t))], which are however 
also not square integrable. Here by "damped free particle" case it is mean 
the dissipative motion under a null external potential. 

This paper is organized as follows. In Section II, we present the Kostin 
model in hydrodynamic form. In Section III, a time-dependent Gaussian An- 
satz for the solution of the Kostin equation is proposed. After expanding the 
arbitrary external potential near the center of mass of the wave packet, it is 
found that the center of mass solves the same Newtonian equation as the one 
for a particle under the external potential and a damping force proportional 
to the velocity. In addition, the width of the wave packet is shown to solve 
the damped Pinney equation [5], [T5], a nonlinear equation whose general 
properties are not fully understood until now. In Section IV, the free particle 
case is analyzed in detail by means of perturbation theory. In Section V the 
obtained solutions are discussed in terms of the Wigner function formalism. 
Section VI is dedicated to the conclusions. 
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2 Kostin model in hydro dynamic form 



In 1972, Kostin [H] proposed the following logarithmic nonlinear Schrodinger 
equation to represent dissipative physical systems, 



<9$ Ti 2 d 2 ^ 

dt 2m dx 2 



His { ^ \ TlV I / * \\ 

V - i — In — + i — ( In — J) 



*, (1) 




where \l/ = ty(x,t) and V = V(x,t) are, respectively, the normalized wave 
function and the external potential, Ti is Planck's constant divided by 2ir, m 
is the mass and u is a positive damping parameter corresponding to Ohmic 
dissipation. For simplicity, in this work only one- dimensional problems are 
considered. Moreover, in Eq. (CQ) we have 

*|, (l , t )|. to (*H), w 

a term introduced in such a manner so that the expectation value for the 
energy becomes just the sum of the kinetic and potential parts. A potential 
term associated to a random force due to interaction with a heat bath could 
have been also added, as in the original Kostin formulation [2]; for simplicity, 
here this possibility is disregarded. 

We introduce the Madelung [TT] decomposition 

* = Aexp(iS/k) , (3) 

where the amplitude A = A(x, t) and the phase S = S(x, t) are real functions, 
and also the probability density p and the quantum fluid velocity v according 
to 

= A 2 v = —— U) 
^ ' m dx 

From this prescription and separating the real and imaginary parts of the 

Kostin equation (CQ) gives 

i+5« - »■ ^ 

dv +y dv 1 dV ^ | h 2 d ( d 2 ^p/dx 2 \ 

dt dx m dx 2m 2 dx \ ^Fp J ' 

which are formally the same as resp. the continuity and force equations 
for a fluid acted by an external potential V, also including a damping term 
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proportional to the velocity and a quantum Bohm potential term (the ~ D 2 
contribution in Eq. (El)). 

In terms of the hydrodynamic formulation, the existence of undamped 
solutions \1/ = <f n (x) exp(—iE n t/H), where <p n (x) and E n are resp. the station- 
ary states and eigenvalues of the Hamiltonian operator H = — (H. 2 /[2m])<9 2 / dx 2 + 
V(x) in the autonomous V = V(x) case, can be understood since these so- 
lutions have zero quantum fluid velocity and hence suffer no friction. 



3 Expansion around classical trajectories 

We propose the following Gaussian Ansatz for the probability density, 

1 ( (x-q) 2 \ m 



where a = a(t) and q = q(t) are auxiliary functions of time to be determined 
in what follows, interpreted resp. as the width and the center of mass of the 
wave packet. Using Eq. ([7]), the only velocity field satisfying the continuity 
equation with a vanishing probability current density pv at infinity turns out 
to be given by 

v = -(x - q) +q . (8) 

Gaussian wave packet solutions to the Kostin equation have been dis- 
cussed by many authors, specially in the harmonic confinement case [2], [0], 
[7], [12]. In this work we highlight the generic applicability of this approach, 
in the following sense. Consider the Taylor expansion 

V(x, t) = V(q(t),t) + V'(q(t),t) (x - q(t)) + _ q{t )f + . . . , ( 9 ) 

where a prime denotes derivative with respect to the first argument. Retain- 
ing up to quadratic terms in V, substituting in the force equation §6§ and 
separating terms proportional to (x — q)° and (x — q) gives 

q + uq = _Iy'(g(t),t), (10) 

Ti 2 

a + ua + u 2 (t)a=- — — , (11) 

4m z a a 
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where we have defined 0J 2 {t) = V"(q(t),t)/m. Hence the center of mass q can 
be identified with the classical particle under the potential V and a linear 
in velocity friction. It should be stressed that Eq. ffTUl) is just the Ehrenfest 
theorem for the potential and distribution (J7|). On the other hand the 
width a satisfy the damped Pinney equation (fTTj) . In the conservative [y = 0) 
case the (nonlinear) Pinney equation is well known to be exactly solvable 
in terms of the solution of the (linear) time-dependent harmonic oscillator 
equation |15j . Approximate solutions exist [5] for weak damping and slowly 
varying non- vanishing frequencies u(t). 

The eikonal S can be found integrating Eq. (J3J) with velocity field given 
by Eq. ([HD, which then gives the wave function 



(2vra 2 ) 1 /4 



exp 



(x — q)' 
4a 2 



exp 



im ( a , 2 .. 



;i2) 



An extra arbitrary function of time could have been added to S, but we 
disregard it since it does not change neither the probability density nor the 
probability current. 

The generic character of the approximate solution (fT2"j) results from the 
arbitrariness of the external potential, as long as the expansion near classical 
trajectories ([9]) holds. The remaining task is to evaluate q, a from Eqs. ( 1101) 
and ( 11 IB . which can be done only in special cases. Curiously, even the free 
particle (V = 0) case does not seems to be solvable in closed form. In the 
next Section a perturbative solution is build for this situation. 



4 The free particle case 

Assuming V = 0, Eq. (llUp is trivially solved as 

<? = Qo + — [1 - exp(-ut)} , (13) 

where q(0) = q , q(0) = %. On the other hand, the solution for the width 
equation 

' a+vh = ^k? (14) 

apparently can not be analytically found. Moreover it can be checked that 
the simple approximate expressions from Ref. [5] developed for a non-zero 
harmonic confinement did not apply in the free case. 
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The ultimate reason why Eq. (j!4p is hard to solve, even approximately 
for small damping, is because the unperturbed solution (for v = 0) grows 
without bound as t — > oo. Indeed, in this conservative case the general 
solution is 



a 2 = (a + d t) 2 + '" " 2 , (v = 0) (15) 



2+2 



ret 



Am? 



where a(0) = a , d(0) = d . 

On the other hand, in the damping dominated case where we can neglect 
the acceleration in Eq. (I14|) and integrating the resulting first-order ordinary 
differential equation, we get 



m?v 



(16) 



for large times, which is consistent with results found from nonlinear theories 
of quantum Brownian motion [IB], [TS]- As pointed out in Ref. [H] in the 
analysis of wave packet spreading in a viscous medium, the leading edge 
experiences a larger friction force than the interior of the Gaussian. In this 
way a retardation of the spreading takes place, which explains the different 
power laws a ~ t for the non-dissipative case in Eq. (|15p and a ~ t 1 / 4 for the 
dissipative case in Eq. f[T5"|) . Here it is proposed a different framework for 
this point, in terms of the collision-dominated regime of the damped Pinney 
equation without harmonic term. 

On general grounds, without neglecting the acceleration term, we have 




= -2vd 2 <0. (17) 

Since every term inside brackets is positive definite, the kinetic energy term 
as well as the inverse square term tend to be negligible for large times. This 
is consistent with the asymptotic result (fT6"j) . Moreover, from Eq. ([TBI) we 
find 

, (18) 

va Aut 

further justifying the neglect of the acceleration term as long as vt ^> 1. 
To start developing a perturbation solution, first rescale according to 

t->z/t, a -> ( — J a, (19) 
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so that Eq. (j!4p is rewritten as 



ed + d = ^- 3 , (20) 

where e is a perturbation parameter eventually set to unity at the end of 
the calculation. The same symbols are used for the original and dimension- 
less quantities, for simplicity of notation. The parameter e is introduced in 
Eq. (T2Uj) . in accordance with our previous considerations indicating that the 
acceleration term can be viewed as a perturbation, except for a transient 
time. 

We then define 

a = a + eai + 0(e 2 ) , (21) 

where a , ai are functions to be determined. Inserting in Eq. f[20"j) and 
equating terms of equal power of e we get 

= - 5< " (23) 

For simplicity we restrict to first order. Higher order corrections are easily 
found, if necessary. 

Equation (1221) is a first order separable ordinary equation, with general 
solution 

a» = (t + c\) l l\ (24) 

where C\ is the integration constant, left indeterminate for the moment. 

Given qj , Eq. (1251) is a first order linear inhomogeneous equation for a\. 
Using Eq. (1241) . the solution is 

c 2 + (3/16) ln(t + cf) 

ai = — ^!F^ — ' (25) 

where c 2 is a further integration constant. 

In conclusion, setting e = 1 we find the approximate general solution 

a = a „ +ai = (t+e ;r + fi±|W±f!) . (26) 
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«(°) = 73-^y + T7-7' (28) 



The two numerical constants c\ and c 2 are determined from the initial con- 
ditions, 

(n\ . ° 2 , 31nCl -fri\ 3 - 12c 2 + 4cf - 91nci 
a(0)=c 1 + - ! + ^- r , a(0) = ^ , (27) 

where we assumed C\ > 0. It turns out that the system (12 7p can be numeri- 
cally solved only. 

Equation f[2"6T) confirm the asymptotic behavior a ~ t 1//4 for large £. Fur- 
ther corrections can be constructed extending the method to second order 
but the result is increasingly awkward without significant gain. 

Eliminating c 2 from Eq. (12 7p results in 

1 3a(0) 3 
4 ~ 4c| 16c} 

which may admit several real solutions after specifying a(0),d(0). For in- 
stance, when a(0) = 2, a(0) = we get C\ = 1.437 and then c 2 = 1.399. The 
numerical solution of Eq. f[2"U|) for these values can be checked to be virtually 
indistinguishable from the perturbation solution (I26|) . However, if the sec- 
ond admissible set c\ = 0.591, c 2 = 0.685 is used, the first-order perturbation 
solution gives not the same excellent agreement with the numerical solution, 
although reproducing the asymptotic behavior. 

Notice that by construction our approximate solution is adapted to the 
collision dominated regime, with the acceleration term treated as a pertur- 
bation. Hence for certain initial conditions one may have some transient 
disagreement between the proposed form and the numerical or the actual 
exact solutions. This is specially true for an initially rapidly expanding wave 
packet, where the effects of the damping become dominant only after a rel- 
atively long time scale. These considerations are verified through the direct 
numerical solution of the damped Pinney equation, with or without harmonic 
confinement term 



5 Wigner function 

It is interesting to briefly consider the Wigner function [23] , which is a valu- 
able tool for the direct interpretation of the quantum system, on phase space. 
For a pure state the Wigner function / = f(x,p,t) is defined as 
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where p is the momentum variable. From Eq. (|2T?|) we have / dxdpf(x, p, t) = 
1. Using the Gaussian form in Eq. (1121) we derive 



/ = -r exp 



p — mq (x — q)\ 



(x — q) 2 2a 2 ( . ma 1 

~la 2 1? 



(30) 



which is valid near classical trajectories for any external potential as long as 
the pair q, a satisfy Eqs. QUI) and (TTT1) . Notice that this Wigner function is 
positive definite and has a Gaussian shape as well. 

From Eq. (|3"0"1) it is tempting to adopt the canonical transformation 

ma 

x — > x — q , p p — mq (x — q) , (31) 

so that in the transformed variables the level curves of the Wigner function 
become ellipses with semi-axis proportional to a and Ti/a. Hence the associ- 
ated area is constant and proportional to h, irrespective of the time evolution 
of the width a. However, clearly these results are a result from the specific 
class of the wave functions, independently of the dynamics generated by the 
external potential. 

In passing, we note that the standard definitions (Ax) 2 =< x 2 > — < 
x > 2 and (Ap) 2 =< p 2 > — < p > 2 can be applied, with averages defined as 
in Eq. ([2]). Restoring physical coordinates the result is 

h 2 

(Ax) 2 = a 2 , (Ap) 2 = m 2 a 2 + — , (32) 

which is in accordance with the uncertainty principle since AxAp > h/2. 
This is an advantage over other approaches using a linear, time-dependent 
Hamiltonian to describe dissipation, like the Kanai-Caldirola model [3], [S], 
which is well-known to violate the uncertainty principle pQ. 



6 Conclusion 

Assuming a Gaussian solution for the Kostin equation and expanding near 
classical trajectories, in this work we have shown that the center of mass of 
the wave packet follows the classical damped dynamics under the external 
potential. Similarly, the width of the Gaussian is shown to satisfy Eq. (II ip . 
which is a non-conservative Pinney equation. Although closed form solutions 
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can perhaps exist for special classes of external potentials, there are well- 
known obstacles to the integrability of the damped Pinney equation: it does 
not match the Tresse-Cartan conditions and hence is not linearizable under 
general point transformations; it does not possess quadratic in velocity cons- 
tants of motion; it can be reduced to an Emden-Fowler equation of index —3, 
but with time-dependent coefficients so that it does not match known solvable 
cases; the use of generalized Sundman transformations may also be proven to 
be useless. See [5] for a more detailed account on the non-solvability of Eq. 
( TTTj) using diverse techniques. Hence in Section IV a perturbation approach 
in the friction dominated case is used, in order to derive for the first time 
analytical approximate solutions of the Kostin equation under zero external 
potential. We note that the same approximate expressions apply in the case 
of a purely time-dependent external force F(t), because in this situation one 
has V(q(t),t) = -F(t)q(t), u(t) = so that Eq. ([II]) still holds. To conclude, 
the general structure of the solutions of Eqs. (1T0|) and ffTTl) . which are gene- 
rically representative in the Kostin model near classical trajectories, remains 
presently largely unexplored. 
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